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Optical isolators will play a critical role in next-generation photonic circuits, but their on-chip integration 
requires miniaturization with suitable nonreciprocal photonic materials. Here, we theoretically 
demonstrate the thinnest possible and polarization-selective nonreciprocal isolation for circularly polarized 
waves by using graphene monolayer under an external magnetic field. The underlying mechanism is that 
graphene electron velocity can be largely different for the incident wave propagating in opposite directions 
at cyclotron frequency, making graphene highly conductive and reflective in one propagation direction 
while transparent in the opposite propagation direction under an external magnetic field. When some 
practical loss is introduced, nonreciprocal isolation with graphene monolayer still possesses good 
performance in a broad bandwidth. Our work shows the first study on the extreme limit of thickness for 
optical isolation and provides theoretical guidance in future practical applications. 

Optical isolators, the photonic analogue of diodes, allow nonreciprocal unidirectional light transmis- 
sion^"^. They have become indispensable components in optical telecommunications and laser applica- 
tions for stabilizing propagating photons and preventing undesired feedback^"^. However, the main 
component of the optical isolators, Faraday rotator, is usually composed of magneto -optical materials with weak 
Faraday effect. It requires a long propagation distance in the order of operating wavelength to achieve good 
isolation, and therefore, optical isolators are generally bulky \ In attempt to reduce costs and realize on-chip 
integration, there is a keen interest to miniaturize the nonreciprocal devices through designing novel structures^. 
Recent advance in fabrication of ultrathin two-dimensional (2D) graphene down to monolayer may provide the 
possibility to further push the optical isolation to atomic scale and has attracted more and more people's interests. 
This is due to the fact that graphene can become nonreciprocal and gyrotropic under an external magnetic 
field^"^^. However, because of its ultrathin thickness, graphene is almost transparent for light^^'^^, and can only 
turn the polarization of incident light by several degrees under a modest magnetic field^^'^^. Hence, external bulky 
resonating structure and polarizer are still needed when incorporating graphene into the traditional isolators^^'^^. 
Achieving optical isolation in the atomic scale is still very challenging. 

In this paper, we theoretically proposed the idea of atomically thin nonreciprocal optical isolation for circularly 
polarized waves by using graphene monolayer. We show that under an external magnetic field, electrons in 
graphene would experience distinct magnitudes of Lorentz force for left-handed or right-handed circularly 
polarized wave (LCP or RCP) incidence at cyclotron frequency. Therefore, for the same circularly polarized 
wave, graphene can be tailored to be transparent to the forward propagation direction, whereas to be reflective to 
the backward propagation direction. Moreover, for the same propagation direction, graphene can be tailored to be 
transparent to one of the circularly polarized waves, whereas to be reflective to the other polarization handedness. 
As a result, the polarization-selective isolator can isolate and filter circular polarizations simultaneously. Our 
work provides an analytical study on the extreme limit of thickness for optical isolation, which might be useful in 
the quest to miniaturize nonreciprocal photonics in the future. 

Results 

Fig. 1 schematically shows the atomically thin and polarization- selective optical isolation based on graphene 
monolayer. Under an external magnetic field perpendicular to graphene plane, graphene (in the xy plane) 

_ / e % 0\ 

becomes gyrotropic and can be described by an asymmetric permittivity tensor: 8= I —iSg 8 0 . For 

\ 0 0 8, 
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Figure 1 | Schematics of the atomically thin and polarization- selective optical isolation for circularly polarized waves. An external magnetic field is 
applied perpendicular to graphene plane. Only one circularly polarized wave can propagate through the gyrotropic graphene in one propagation 
direction. The arrows in the helixes show the instantaneous electric field distributions of circularly polarized waves. 



plane wave propagating along magnetization {±z) direction, the 
eigenmodes in gyrotropic graphene are LCP and RCP. Without loss 
of generality, we analytically calculate the transmissivity of circularly 
polarized waves normally incident from a dielectric medium with 
relative permittivity £i onto the gyrotropic graphene layer, and trans- 
mitted into a second dielectric medium with relative permittivity £3. 
Detailed calculation can be found in Methods. The calculated wave- 
vectors for LCP and RCP in gyrotropic graphene are as follows: 



j^+RCP = co/c^8 — 8g for forward propagation ( + 

j^-RCP = (jrf/cy^8-\-8g for backward propagation ( — z), 

from which we see that for the same circularly polarized wave they 
have different phase velocities {Vph = co/k) in the forward and back- 
ward propagations. Moreover, since |A:"^^^^| = lA:"^*^^! and 
I ^ -LCP I ^ |^+i?cP|^ ^YiQ transmissivity of LCP in one propagation dir- 
ection is the same with that of RCP in the opposite propagation 
direction. With this symmetry, we only need to analyze the forward 
propagating transmissivity. 

The ideal case of graphene monolayer without loss is firstly dis- 
cussed in Fig. 2(a). The parameters in Fig. 2(a) are as follows: the 
chemical potential fi^ = 0.5 eV, the cyclotron frequency wJ2n = 
10 THz, the damping constant y = 0, and Si = £3 = L From the 
analytical derivations (details are in Methods), we can get forward 
propagating transmissivity t^^^^ = 0 and t^^^^ = 0.99 at cyclotron 
frequency, indicating good isolation performance. The underlying 
mechanism of the remarkable contrast in transmissivity at cyclotron 
frequency is that the electron velocity in graphene under LCP incid- 
ence is largely different from that under RCP incidence, making 
gyrotropic graphene highly conductive and reflective for one circular 
polarization incidence while transparent for another circular polar- 
ization incidence. From a microscopic point of view, by using free 
electron modeP^, niejfdv/dt = — e(E-\-v x 5), where m^jf is the effec- 
tive electron mass and E is the electric field of incident wave being 
either RCP or LCP. Under time harmonic excitation, we can derive 



- e [icoE -{-cOcXE codcoc 'E)] 

CO 



(2) 



where cDc= —eB/nieff. In this paper, we set the external magnetic 
field B along — z direction, and the scalar cyclotron angular frequency 



is defined as C0c=\cbc\. For the forward ( + z) propagation, we have 

-ieE^^^^ 

for LCP incidence with E^^^^ = Eo(x-iy) 



VLCP _ 



and 



mejf(co-Wc) 

_:^p + RCP 

RCP _ 



for RCP incidence with E^^ 



meff(w + Wc) 

Eo(x-\-iy). Hence, the magnitude of Lorentz force acting on graphene 
electrons under LCP incidence (|F^^^^| = |— — 1) is much 



+ RCP 



(W-Wc) 

-ewE+^^\ 



at 



larger than that under RCP incidence (|F , , . 

CO = (Dc. As a result, vtjf^^ can be theoretically very large and vt^^^^ is 
reduced to half at cyclotron frequency. Since the electric current 
density Jocv and J = aE, gyrotropic graphene becomes highly con- 
ductive under LCP incidence, indicating gyrotropic graphene behaves 
like a perfect electric conductor, and LCP is thus totally reflected. In 
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Figure 2 | Transmissivity of RCP (f-'^^O and LCP (f-'^^O through 
gyrotropic graphene monolayer in the forward propagation, (a) 

Comparison of f^^^^ and t^^^^ under the lossless {InylcOc = 0) and loss 
{Iny/Wc =1) conditions with wJ2n =10 THz. (b) Loss influence on the 
isolation performance with different cyclotron frequencies. f-^^cpi^+LCP 
increases as 2ny/Wc or co^ decreases. Other parameters are as follows: lu^ = 
0.5 eV and 81 = 83 = 1. 
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contrast, for RCP incidence, graphene has a finite negative permittiv- 
ity and wave is evanescent inside of the graphene. As its skin depth 
(l/Im(A;^^^0 > 100 nm) is far larger than graphene's atomic thick- 
ness, graphene is highly transparent to RCP incident wave. 

Because of the collision between electrons, loss is inevitable 
in practical graphene samples. In the following part, we studied 
the loss influence on the isolation performance. When loss is 
considered, the electron velocity in equation (2) becomes 

— e[/(ft) + iy)E + X £ — ^ 

V = 2 ~ • 1^0^ the forward ( + z) 

nieffiico + iy) -w^^) 

propagation and by setting B along 

meff(co + iy-(Dc) 

; for RCP incidence. Therefore, as long as InyloJc 

mejf{co + iy + cOc) 

is sufficiently small, v^J^^^ can still be theoretically very large, and 
VB"^o^^>VB"t?^^ at cyclotron frequency can still be achieved. Fig. 2(a) 
shows the analytical transmissivity of LCP and RCP under the loss- 
less {Inylojc = 0) and loss {Inylojc = 1) conditions. One can see that 
the loss degrades the performance of isolation. Fig. 2(b) shows the 
transmissivity ratio between RCP and LCP at cyclotron frequency 
(j.+RCPi^+LCP^ as a function of InylcJc with different cyclotron frequen- 
cies, from which one can see that, on one hand, t^^cpi^+icp incj-g^ses 
as Inylojc decreases, on the other hand, with the same value of Inyl 
increases as the cyclotron frequency decreases. Hence, 



RCP I ^+ LCP 



better isolation performance can be obtained by keeping InylcOc 
sufficiently small or setting smaller cyclotron frequency. 

To verify the analytical prediction of the atomically thin isolation. 
Fig. 3 shows the numerical simulations of the lossless case at 10 THz 
with the use of Finite Element Method (COMSOL Multiphysics). In 
the simulation, the following parameters are used: — 0.5 eV, ojjln 
= 10 THz, y = 0, and Si = £3 = 1. We set the incident field with unit 
magnitude. When propagating in the forward direction, RCP is 
almost totally transmitted through graphene, as shown in Fig. 3(a), 
featured with unit magnitude of maximum Ey at both sides of gra- 
phene. Meanwhile, LCP is totally reflected, as shown in Fig. 3(b), 
characterized with twice unit magnitude of maximum Ey at the left 
graphene side and near zero at the right side. Moreover, from 
Fig. 3(b), one can see that the total electric field at the graphene 
boundary is always zero, indicating that graphene exhibits zero 
impedance, verifying our previous theoretical analysis that in this 
case graphene behaves like a highly conductive material. When prop- 
agating in the backward direction, RCP is totally reflected in Fig. 3(c) 
and LCP is almost totally transmitted in Fig. 3(d), behaving inversely 
to the forward propagation. 
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Figure 3 | Simulation demonstration of the atomically thin optical 
isolation for circularly polarized waves under the lossless condition. Ey 

field distribution of RCP and LCP waves propagating in (a-b) the forward 
( +z) and (c-d) backward ( — z) directions at 10 THz. The black dashed 
lines denote the location of graphene and the green arrows represent the 
wave incident direction. The parameters are: fic = 0.5 eV, ajjln = 
10 THz, y = 0, and £1 = £3 = 1. Under the lossless condition, graphene 
monolayer can isolate and filter the circular polarizations simultaneously. 



The numerical verifications for the loss case are shown in Fig. 4. 
Based on previous practical graphene parameters, we set cojln = 
0.3 THz, fic = 0.5 eV, and y = 0.63 THz in the demonstration. These 
correspond to B = —1 Tesla, the 2D carrier density Us = 2 X 
10^^ cm"^ and the electron mobility ^ = 28000 cm^V^s"^ These 
parameter settings are achievable in experiments, because that the 
2D carrier density Us in graphene is experimentally reported to be 
gated efficiently from very low (10^° cm"^) to very high values 
(10^^ cm~^y^~^^, and the electron mobility up to ^ = 
60000 cm^V^s"^ has been experimentally obtained in high-quality 
suspended^^ or hexagonal boron nitride substrate supported^^ gra- 
phene samples. We also consider the graphene is mounted on a 
substrate. The influence of surrounding dielectric substrate on the 
isolating performance is shown in Fig. 4(a,b), where symmetric 
structure is assumed with £1 = £3 = to avoid additional back- 
ground reflection. Smaller £^ is found favorable to achieve better 
isolation performance. To ensure the high electron mobility in 
graphene and therefore a smaller damping constant y, hexagonal 
boron nitride with = 4.8^^ is chosen. From Fig. 4(b), one can 
see that high t^^cpj^+LCP obtained in a wide frequency range 

with the maxima of 25.96 at 0.311 THz (3.7% higher than co^ = 
0.3 THz). When t^RCP^+icp > jg defined as a reference, a 
0.26 THz bandwidth can be obtained with frequency ranging 
from 0.2 to 0.46 THz. Fig. 4(c-f) show the numerical simulations 
of the nonreciprocal isolation at 0.3 THz. In the forward propaga- 
tion, RCP wave passes graphene with a transmissivity of 31.5% in 
Fig. 4(c), while LCP is efficiently blocked with a trivial transmissivity 
of 1.23% in Fig. 4(d). When propagating in the backward direction, 
it behaves inversely in Fig. 4(e,f). From these demonstration, one see 
that an atomically thin wide-band THz nonreciprocal optical isola- 
tion is possible by using practical graphene samples, showing the 
appealing potential applications of graphene in ultrathin nonreci- 
procal device design. 
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Figure 4 | Atomically thin optical isolation under the loss condition, (a) 

Infiuence of on the forward propagating transmissivity of LCP and RCP 
waves, where Si = £3 = (b) The isolation performance represented by 
^+RCPi^+LCPf^^ different e^. (c-f) Enfield distribution of circularly polarized 
waves propagating in the forward ( + z) and backward ( — z) directions at 
0.3 THz. The Parameters are: fic = 0.5 eV, C0c/2n = 0.3 THz, y = 
0.63 THz, and = 4.8. The black dashed lines denote the location of 
graphene, and the green arrows represent the wave incident direction. 
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Discussion 

Previous graphene based traditional Faraday isolators^^'^^ are gen- 
erally composed by a graphene sheet, a resonating structure and two 
wire grids. In that Faraday isolator structure, the gyro tropic graphene 
only serves as a Faraday rotator, the thickness of the resonating 
structure is around half of the operating wavelength^^, and the two 
wire grids serving as the polarizers are indispensable to realize isola- 
tion. When considering the thickness of the resonating structure and 
wire grids, the total thickness of previous graphene based Faraday 
isolators^^'^^ is usually in the order of operating wavelength, far larger 
than our atomic thickness. Hence, previous graphene based Faraday 
isolators^^'^^ are different from the demonstration in this work where 
no external wire grids or other resonant structures are needed. 
Flexible tunability is a potential advantage of graphene based 
devices^^'^^. Since the best performance of the ultrathin graphene 
isolation exists near cyclotron frequency, the isolation working fre- 
quency range should also be tunable by external magnetic field. 

In conclusion, we have theoretically proposed the atomically thin 
and polarization -selective nonreciprocal optical isolation for circular 
polarizations. The mechanism underlying the atomically thin optical 
isolator is revealed from a microscopic viewpoint. Our work shows a 
feasible solution to control the propagation of circularly polarized 
waves in graphene and to manipulate light nonreciprocally at atomic 
scale, which opens up new possibilities for further innovations in 
functional ultrathin nonreciprocal optics. 

Methods 

Graphene monolayer is treated as a homogeneous film with thickness d = I nm^°'^^ 
and characterized by an anisotropic permittivity tensor diag[£e<j Seq ej . Based on Kubo 
formula, Sea can be cast to the Drude modeF'^"'^"'^*^'^'', 



= 1 + /- 



a + Gin 



CQEod 



COEod 



cq{cq + iy) 



(3) 



where Gin 



+ 2ln 



-^iJkBT J 



+ 1 ) I and Ginter are optical con- 



ductivity attributed to intra-band and inter-band transitions respectively, Ej^, 



iGinter^CO^od, (Dp 



+ 2ln 



1/2 



is the plasma fre- 



jiH^Eod yksT 

quency, y is the inverse of relaxation time, and £o> ^> ^> and /c^ are the dielectric 
permittivity in vacuum, electron charge, reduced Planck constant, and Boltzmann 
constant, respectively. When setting fic = 0.5 eV and temperature T = 300 K, we 



have the effective electron mass f 



-- 0.1 X 0.91 X 10-'° kg and (Oplln = 410 THz. 



At the operation frequency below 1 THz, Einter is equal to 3.9. 

Gyrotropic graphene can be described by an asymmetric permittivity tensor 

E iEg 0 \ 

— iEg E 0 I , where e = Sinter — (c^ + ^ V) / ['^(('^ + ^ 7 )^ — c^c )] 
0 0 eJ 

Eg= —cOpCOc/ [co((co + iy)^ — c«c)] • Hence, e + £^ = £inter — cOp/{co(co — a>c + iy)) and 
E — Eg = Sinter — (^p/ (co(aj + a»c + iy)). The scalar cyclotron angular frequency cOc is 
defined as |ft)c|, when the external magnetic field B is along —z direction. 

Wave polarizations can be viewed by either the temporal or spatial view points. In 
the case of a circularly polarized wave, from the temporal view point, the tip of the 
electric field vector describes a circle at a fixed point in space as time progresses. If the 
right-hand (left-hand) thumb points in the direction of propagation while the fingers 
point in the direction of the tip motion, the wave is defined as right-hand (left-hand) 
circularly polarized. From the spatial view point, the instantaneous electric field 
vector of the RCP (LCP) wave describes a left-hand (right-hand) helix along the 
direction of propagation, as shown in Fig. 1. 

To achieve the condition of good isolation for circular polarizations, the trans- 
missivity of the forward propagation in Fig. 1 are analytically calculated in the fol- 
lowing. The calculation of backward propagation can be done by the same 



methodology. Setting 



-L/RCP _ 



L/RCP Jk+'^^^'z 



-c 



L/RCP -i/c+^Z^^'^, 



l^'^ region with / = 2 standing for gyrotropic graphene and / = 
dielectric. After calculation, we have: 



zj>) in the 
1,3 for isotropic 



^L/R. 



j^L/RCP i2k;-';'''''d,, ,,f^L/RCP ul/RCP. 



.L/RCP ,-2Lt'/"^" 



'da , rfL/RCP.^L/RCP / .L/RCP^ 



(4) 



where R, 



L/RCP _ 



L/RCP 
"/(/+!) rj^L/RCP _ 
L/RCP' ^1(1+1) ~ 



1(1+1) 



^_ L/RCP _ 

L/RCP ' 



kl{E±Eg) 



^ L/RCP 



dzi = 0, and dz2 = d. The reflection and transmission coefficients are 
j^+L/RCP^cL/RCP^^L/RCP y+v«cp^^V«cP/^t/fiCP_ ^^.p^.tively. After some 
algebra, we have 



^ L/RCP _ 



Sinter 



inter 



oj(oj -\-0Jc-\- iy) 



+ + iy) 



Sinter 



Sinter 



oj{oj -f a>c + iy) 



+ a»c + iy) 



(3) 



1 + - 



\/sr— \ / Sinter 



(i)(co -I- a>c + iy) 



\/sr+ A /Sinter 



+ a>c + iy) 



a){a)-\-a)c-\-iy) ^ak^'^^'^'d 



+ a»c + iy) 



J.L/RCP ^L/RCP Uk^^l^'^'' -k,)d 
j^ + L/RCP _ ^12 ^23 ^ 

h/K 



2V^ 



\/Sr+'\/ Sinter 



Sinter 



Sinter 



-I- a>c + ^y) V "^'^^^'^ a>(a> -I- + zy) i2k^^^^'^^d 



(6) 



\ I Sinter } — =r-^ — \ I Sinter } — =r-^ — 

09(09 + CL»c + zy) y cl»(cl» + o^c + zy) 

The reflectivity and transmissivity are r+^/^^^ = |K+^/^cp|2 
^+L/RCP = |7^+V-RCP|2^g^^g^^ respectively. From equation (5) and equation (6), 
^+LiRCP ^+LiRCP Strongly dependent on oOc, cOp, and y, when the operating 
frequency is around cOc. In particular, for LCP incidence in the lossless case, when co 



, |£ + e^| = |£inter-COp/(co(cO-C«c))H0O, |pff^HoO, \p\^^\-- 



^ . =0,i^ff^=-l, 

j^LCP ^ ^ jLCP ^ Q ^LCP ^ 2, and thus r+^^^ = 1 and ^^^^ = 0. For RCP incidence in 

23 12 23 

the lossless case at co = cOc, E — Eg = Sinter — l{oj{(D + (Dc)) is negatively large but 
finite, and we get t^^'^^ = 0.99 as shown in Fig. 2(a) through numerical computation. 
When considering the loss case and setting oOp as a constant, t^^'^^^ is strongly 
dependent on znylcOc and cOc, which can be seen from Fig. 2(b). 

When /ic and y are given, the corresponding Us and /i can be theoretically calculated 
by using ns = pi^J (nH^ v^Y'^'^'^ and l/y = ^^^/{ev^y^'^'^, respectively, where Fermi 
velocity Vp ~ 0.95 X 10^ m/s. When \B\ = 1 Tesla and fic = 0.5 eV, we can neglect 
Landau quantization influence and use the Drude model directly for graphene 
according to Ref. [10], because the following conditions are satisfied: the Landau 

energy scale L = (iHeBvl) = 0.035ev <C fic and the operating frequency (set below 
1 THz) energy HcQ<^fi^. 
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